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(57) Abstract 

A receiving apparatus Doppler spread estimator includes a receiver receiving signals over a transmission channel that varies over time 
due to fading and converting the signals to digital samples for processing. A stationarizer operatively associated with the receiver stationarizes 
the digital samples. An estimator is operatively associated with the stationarizer for estimating spectral density of the stationarized digital 
samples. The correlator is operatively coupled to the estimator for correlating the spectral density estimates with functions corresponding 
to hypotheses of plural Doppler spread values to select one of the Doppler spread values as an estimate of actual Doppler spread. 



. <WO 005756SA 1 _l_> 



FOR THE PURPOSES OF INFORMATION ONLY 



Codes used to identify States party to the PCT on the front pages of pamphlets publishing international applications under the PCT. 



AL 


Albania 


ES 


Spain 


LS 


Lesotho 


SI 


Slovenia 


AM 


Armenia 


FI 


Finland 


LT 


Lithuania 


SK 


Slovakia 


AT 


Austria 


FR 


France 


LU 


Luxembourg 


SN 


Senegal 


AU 


Australia 


GA 


Gabon 


LV / 


Latvia 


SZ 


Swaziland 


AZ 


Azerbaijan 


GB 


United Kingdom 


MC 


Monaco 


TD 


Chad 


BA 


Bosnia and Herzegovina 


GE 


Georgia 


MD 


Republic of Moldova 


TG 


Togo 


BB 


Barbados 


Gil 


Ghana 


MG 


Madagascar 


TJ 


Tajikistan 


BE 


Belgium 


GN 


Guinea 


MK 


The former Yugoslav 


TM 


Turkmenistan 


BF 


Burkina Faso 


GR 


Greece 




Republic of Macedonia 


TR 


Turkey 


BG 


Bulgaria 


HU 


Hungary 


ML 


Mali 


TT 


Trinidad and Tobago 


BJ 


Benin 


IE 


Ireland 


MN 


Mongolia 


UA 


Ukraine 


BR 


Brazil 


IL 


Israel 


MR 


Mauritania 


UG 


Uganda 


BY 


Belarus 


IS 


Iceland 


MW 


Malawi 


US 


United States of America 


CA 


Canada 


IT 


Italy 


MX 


Mexico 


uz 


Uzbekistan 


CF 


Central African Republic 


JP 


Japan 


NE 


Niger 


VN 


Viet Nam 


CG 


Congo 


KE 


Kenya 


NL 


Netherlands 


YU 


Yugoslavia 


CH 


Switzerland 


KG 


Kyrgyzstan 


NO 


Norway 


ZW 


Zimbabwe 


CI 


Cdte d'Tvoire 


KP 


Democratic People's 


NZ 


New Zealand 






CM 


Cameroon 




Republic of Korea 


PL 


Poland 






CN 


China 


KR 


Republic of Korea 


PT 


Portugal 






CU 


Cuba 


KZ 


Kazakstan 


RO 


Romania 






CZ 


Czech Republic 


LC 


Saint Lucia 


RU 


Russian Federation 






DE 


Germany 


LI 


Liechtenstein 


SD 


Sudan 






DK 


Denmark 


LK 


Sri Lanka 


SE 


Sweden 






EE 


Estonia 


LR 


Liberia 


SG 


Singapore 







BNSDOCID: <WO 0057568A1J_> 



WO 00/57568 



1 



PCT/US00/05129 



DOPPLER SPREAD ESTIMATION SYSTEM 

FIELD OF THE INVENTION 
This invention relates to wireless communication systems and, more 

particularly, to a method and apparatus for Doppler spread estimation in a wireless 

5 communication system. 

BACKGROUND OF THE INVENTION 
In wireless mobile communication systems the transmission channel can vary 
with time. This is due to the mobility of the transmitter and/or the receiver. Variations in the 
transmission channel are characterized by the Doppler spectrum. The Doppler spectrum is the 

10 Fourier transform of the autocorrelation function of the channel's time impulse response. The 

width of the Doppler spectrum is referred to as the Doppler spread in "Microwave Mobile 
Communications," W. C. Jakes, Jr., ed., New York: Wiley, 1974, IEEE Press, 1994, and 
"Wireless Communication," Theodore S. Rappaport, Prentice Hall PTR, 1996. The Doppler 
spread is commonly used as a measure of the channel variation rate or fading rate. 

15 It is often desirable to estimate the Doppler spread. For example, in order to 

adaptively optimize a coherent detector in a receiver, the parameters of the adaptation 
algorithm are made to be a function of the estimated Doppler spread. Also, in a cellular mobile 
telephone system, the handoff process can be enhanced if an estimate of the Doppler spread 
is available. This avoids handoff of fast moving users to micro cells. 

20 U.S. Patent Nos. 4,723,303 and 5,016,017 discuss the use of Doppler spread 

estimation. Both patents describe estimating the Doppler spread from the received signal itself 
However, these methods may not give a good estimate with a non-stationary received signal; 
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Lars Lindbom, "Adaptive Equalization for Fading Mobile Radio Channels," Techn.Lic.Thesis 
No. UPTEC 92124R, November 1992, Department of Technology, Uppsala University, 
Uppsala, Sweden, discloses a method of estimating the Doppler spread from a sequence of 
channel estimates. This method uses differentials of the channel estimate, which comprise 
5 difference of values between two points in time. However, the differential itself is generally 

very noisy and needs averaging. As a result, the average gives biased estimates of the Doppler 
spread. 

Still another proposal improves Doppler estimation by adaptive filters 
controlled by the final Doppler spread estimate. However, this is done with a significant 
1 0 increase of receiver complexity, while the estimate remains biased. 

The present invention is directed to overcoming one or more of the problems 
discussed above in a novel and simple manner. 



SUMMARY OF THE INVENTION 
In accordance with the invention there is provided a receiver Doppler spread 
15 estimator which utilizes stationarization of a non-stationary received signal. 

Broadly, in accordance with one aspect of the invention, there is disclosed a 
Doppler spread estimator including receive means for receiving signals over a transmission 
channel that varies over time due to fading, and converting the signals to digital samples for 
processing. First processing means are operatively associated with the receive means for 
20 stationarizing the digital samples. Second processing means are operatively associated with 

the first processing means for estimating spectral densities of the stationarized digital samples. 
Third processing means are operatively coupled to the second processing means for correlating 
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the spectral density estimates with functions corresponding to hypotheses of plural Doppler 
spread values to select one of the Doppler spread values as an estimate of actual Doppler 
spread. 

It is a feature of the invention that the receive means is adapted to receive code 
5 division multiple access (CDMA) signals. 

It is a further feature of the invention that channel coefficients are obtained 
directly from the CDMA signals. 

It is another feature of the invention that the receive means is adapted to receive 
signals including pilot symbols spread with a known CDMA spreading code. 
10 It is another feature of the invention that the received signal represents a 

sequence of transmitted symbols and the first processing means comprises a mixer for mixing 
the digital samples with a complex conjugate of the transmitted symbols. The complex 
conjugate of the transmitted symbols is derived from a channel estimator. 

It is another feature of the invention that the first processing means comprises 
15 a low pass filter. 

It is another feature of the invention that the second processing means develops 
vector data representing spectral components. 

It is still another feature of the invention that the third processing means 
comprises a multi -channel correlator. The correlator calculates correlations between estimated 
20 spectral densities and weighting functions corresponding to the hypotheses of the plural 

Doppler spread values. The third processing means includes a selector for selecting the 
hypothesis producing a highest correlation value. 
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In accordance with another aspect of the invention, there is disclosed a 
receiving apparatus Doppler spread estimator including a receiver receiving signals over a 
transmission channel that varies over time due to fading and converting the signals to digital 
samples for processing. A stationarizer operatively associated with the receiver stationarizes 

5 the digital samples. An estimator is operatively associated with the stationarizer for estimating 

spectral density of the stationarized digital samples. The correlator is operatively coupled to 
the estimator for correlating the spectral density estimates with functions corresponding to 
hypotheses of plural Doppler spread values to select one of the Doppler spread values as an 
estimate of actual Doppler spread. 

1 0 There is disclosed in accordance with a further aspect of the invention a method 

of estimating Doppler spread of a received signal comprising the steps of receiving the signals 
over a transmission channel that varies over time due to fading, and converting the signals to 
digital samples for processing; stationarizing the digital samples; estimating spectral density of 
the stationarized digital samples; and correlating the spectral density estimates with functions 

15 corresponding to hypotheses of plural Doppler spread values to select one of the Doppler 

spread values as an estimate of actual Doppler spread. 

More particularly, the present invention relates to a method and apparatus for 
Doppler spread estimation in a wireless communication system. The Doppler spread 
information can be used to optimize performance of a communication system. Applications 

20 include both base and mobile station receivers where a coherent detector is the preferred 

detector type. 

Further features and advantages of the invention will be readily apparent from 
the specification and from the drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
Figs. 1 -3 are block diagrams of receiving apparatus including a Doppler spread 
estimator according to the invention; and 

Fig. 4 is a block diagram of the Doppler spread estimator according to the 

5 invention. 

DETAILED DESCRIPTION OF THE INVENTION 
Referring to Fig. 1, a receiving apparatus 10 according to the invention is 
illustrated. The receiving apparatus 10 is adapted for use with known pilot symbols. The 
receiving apparatus 10 includes an antenna 12 for receiving radio signals represented at 14. 

10 The antenna 12 is coupled to a radio receiver and converter 16. The radio receiver and 

converter 16 filters amplified signals and converts them to a suitable form for processing, such 
as complex numerical sample values. The signals are in the nature of digital samples for 
processing. The radio receiver and converter 16 is coupled to a channel estimator 18r*The 
channel estimator 18 is connected to a Doppler spread estimator 20 and receives known 

15 symbols from a known symbol block 24. The known symbols comprise pilot or reference 

symbols also included in the received signal used for determining channel estimates. The 
known symbols are stored in a memory or generated locally in a code generator. The channel 
estimator 1 8 correlates the received digital samples with the known symbols to provide channel 
estimate to the Doppler spread estimator 20. The Doppler spread estimator 20 estimates the 

20 Doppler spread which is sent to a signal processing block 22. The signal processing block 22 

processes the sampled signal to extract information. The signal processing block 22 is also 
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connected to the channel estimator 1 8 so that after Doppler estimation the channel estimation 
can be improved. 

Referring to Fig. 2, a block diagram of a receiving apparatus 10' is illustrated. 
The receiving apparatus 10' is generally similar to the receiving apparatus 10 of Fig. 1. For 
5 simplicity, like elements are referenced with like numerals, while modified elements are 

indicated with primed numerals. Particularly, the receiving apparatus 10' differs in replacing 
the known symbols at block 24 with a symbol estimator 24'. This receiving apparatus is used 
in an application where symbols are unknown. The channel estimator 1 8 attempts to estimate 
the symbols, which are then stored in the symbol estimator 24'. The estimated symbols are 

10 subsequently used by the channel estimator 18 for determining channel estimates, as above. 

In the United States code-division multiple access (CDMA) cellular system 
known as IS-95, a transmitter transmits a stream of known symbols known as the pilot code. 
The pilot code is transmitted on the same channel at the same time as other, information 
bearing, symbols using different spreading codes. Fig. 3 illustrates a receiving apparatus 10" 

1 5 used in such a CDMA system. With CDMA systems the channel can be estimated directly to 

be used for Doppler spread. In this case, the known symbol block 24 of Fig. 1 is omitted. The 
channel estimator 1 8 correlates the received signal, comprising the pilot code and other codes 
in additive superimposition, with the known pilot code, and low pass filters the resultant 
complex correlation to obtain channel estimates. The received signal is also correlated with 

20 other codes carrying information it is desired to decode. The results of correlating with 

information carrying codes are multiplied by the conjugate of pilot code correlations for the 
same delay, and the results added to coherently combine the multi-path signals. 
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In wide band CDMA (WBCDMA) systems, modulation symbol intervals are 
much shorter. This allows multiple propagation paths to be resolved with much finer time 
resolution. 

In mobile communication systems, the received signal in a flat Rayleigh channel 
5 can be described by the following baseband model: 

r(t) = s(t) ■ c(t) + n(t), 

where 

n 

d n is an n-th transmitted symbol, pf/J is a transmitted pulse, c(t) is a complex Gaussian process 
10 with mean zero and envelope correlation function 

Kc(t) = cr c 2 Jo((Q<*t\ 

&1 is a variance of the process c(t) y J 0 (x) is a Bessel function of the first kind and of order 
zero, 

co d = (OqV/Cq, cd 0 is a carrier frequency, v is a mobile's velocity, c 0 is a light speed, and n(t) is 
1 5 white noise with the spectral density g. 

There are two conventional approaches for the basic Doppler spread estimation 
solution. These approaches are spectral analysis of the received signal or correlation analysis 
of the received signal. With spectral analysis of the received signal: 



i o-i 

1 x — \ 



L nexp(-j2x—) 

i=qN ™ 
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where 

n = di-a + 

r„ c„ and w, are the i-th samples of the processes r(t) 9 c(t) y and n(t) y respectively. For this 
algorithm an estimate of the Doppler spread can be expressed in terms of FFT coefficients: 

r 

5 /rf= -77arg[max5(A:)] 

N 

where f q is a sampling frequency. 

With correlation analysis of the received signal: 

| HQ 

An estimate of the Doppler spread can be expressed in terms of the correlation function 
10 K(m), calculated for some value m 0 : 

where K c "'(x) is a function which is the inverse of the function K c (x). 

However, a vector r = (r,, r N ) is a time-varying, non-stationary random 
discrete process. Traditional spectral analysis does not work for non-stationary processes. To 
1 5 solve this problem, the present invention takes into consideration that d s = exp(j<f>,) Then, the 

algorithm 
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where d i is an estimate of the symbol d i , transforms the non-stationary vector r to the 
stationary vector f = . 

Indeed, for non-dispersive channels di ~ di . Besides, statistical characteristics 



of the noises /*, and hi = di ♦ rii are the same, because 



di 



= 1 . Therefore, 



rt * a + A7i 

It is possible to improve the algorithm (9) by using the moving average: 

j iM+N0-\ 

r<=— Z (d,)V», 



In this case the errors Adn —d — dn are reduced and Equation (10) becomes more precise. 
For CDMA systems, channel coefficients (more exactly, channel coefficients estimates c i ) are 
10 directly obtained. In this case, channel coefficient estimates can be processed instead of r i . 

In accordance with the invention the non-stationary vector r is transformed to 
reduce the problem to analysis of the stationary vector f . Thereafter, the invention uses the 
optimal approach for estimating Doppler spread. This approach uses a criterion of the 
maximum of likelihood ratio. The optimal estimation algorithm is obtained by using the 
1 5 following equation: 
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= 0, 



where A(r\ fd) is the likelihood ratio which is formed based on the vector r . 

Because the vector f has a Gaussian distribution, with the accuracy up to an 
inessential constant, the following equation can be written for the log-likelihood ratio: 



A(f 



«/2-l 



k=0 



where Gff(k\fd) is a weight function which corresponds to some value of the Doppler 

spread f d . 

The method of estimating Doppler spread in accordance with the invention 
utilizes stationarization of the received signal, spectral analysis of the stationary vector, and 
10 correlation of the spectral estimate. Particularly, for stationarization of the received signal: 



No 



Spectral analysis of the stationary vector f = (fi,..., W) using a fast Fourier 
transform (FFT): 
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X fi exp(-;2,r— )| 2 



Correlation is done in the frequency domain of the spectral estimate S(k) with 
L weight functions Gff(k\fd) which correspond to L values of the Doppler spread f d \ 



A(r 



N/2-1 



*=0 



5 where f d = {fji}) y ... 9 fJL)} is a vector of Doppler spread hypotheses. One of the Doppler 

spread hypothesisyy/w^) is selected if: 

A(r|^(m))> A(f\f d (p)) 

Fig. 4 illustrates a block diagram of the Doppler spread estimator 20 of Fig. 1 . 
Stationarization of the received signal is performed using a multiplier 26 and low pass filter 28. 
10 The multiplier 26 receives the sampled received signals r tt which are multiplied with the 

complex conjugate of the transmitted symbols dn . The radio receiver and converter block 

16, see Fig. 1, is the source of the received signal r„. The source of the complex conjugate of 
the transmitted symbols can be the channel estimator 1 8 or the known symbols 24. The output 
of the multiplier 26 is passed through a low pass filter 28 which reduces statistical errors 

1 5 Adn = dn — dn of the transmitted symbol estimation and out-of-band noise. 
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The output of the filter 28 is coupled to an input of a second processing block 
30. The block 30 estimates the spectral density of the incoming signal h using the algorithm 
of Equation (15). The output, represented in bold lines, indicates that the output data is in the 
form of vector data which are the spectral components S = S(N 1 2)] . 

5 The output of the block 30 is fed to an input of a third processing function 32 

in the form of a multi-channel correlator which produces the likelihood ratio metrics using 
Equation (16). Each channel of the correlator calculates correlations in frequency domain 
between the spectral estimate S(k) and the waiting function Gff(k\id) . The waiting function 

G//(Jfc|f d(m)) for the ro-th channel corresponds to some selective value//w) of the Doppler 

10 spread. Particularly, these comprise hypotheses of plural Doppler spread values. The 

likelihood ratio metrics are compared with each other using Equation (16). The outputs are 
applied to a max function block 34. On the basis of comparisons between the likelihood ratio 
metrics of L Doppler spread hypotheses, the multi-channel correlator 32 generates the decision 
that the Doppler spread h//m) if the m-th channel has the maximum output value. 
! 5 The illustrated method for estimating Doppler spread can be used in WCDMA 

or IS-95 systems, since in these systems known symbols or pilot symbols are periodically 
inserted into the data stream at every slot. 

As will be appreciated by one of ordinary skill in the art, the present invention 
may be embodied as methods or devices. Accordingly, the present invention may take the 
20 form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment 

combining hardware and software aspects. The present invention has been described in part 
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with respect to the block diagram illustrations of Figs. 1-4. It will be understood that each 
block of the illustrations, and combinations of blocks, can be implemented by computer 
program instructions. These program instructions, which represent steps, maybe provided to 
a processor to produce a machine. 
5 Accordingly, blocks of the block diagrams support combinations of means for 

performing the specified functions in combinations of steps for performing the specified 
functions. It will be understood that each block of the illustrations, and combinations of 
blocks, can be implemented by special purpose hardware-based systems which perform the 
specified functions or steps, or combinations of special purpose hardware and computer 
1 0 instructions. 

Thus, in accordance with the invention, a Doppler spread estimator utilizes 
stationarization of non-stationary received signals and calculates correlation functions which 
are compared to hypotheses to determine an accurate and reliable estimate of Doppler spread. 
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WE CLAIM: 

1 . A Doppler spread estimator comprising: 

receive means for receiving signals over a transmission channel that varies over 
time due to fading and converting the signals to digital samples for processing; 

first processing means operatively associated with the receive means for 
estimating spectral densities of the digital samples; and 

second processing means operatively coupled to the first processing means for 
correlating the spectral density estimates with functions corresponding to hypotheses of plural 
Doppler spread values to select one of said Doppler spread values as an estimate of actual 
Doppler spread. 

2. The Doppler spread estimator of claim 1 wherein said receive means is 
adapted to receive Code-Division Multiple Access (CDMA) signals. 

3. The Doppler spread estimator of claim 2 wherein channel coefficients 
are directly obtained from the CDMA signals. 

4. The Doppler spread estimator of claim 1 wherein said receive means is 
adapted to receive signals including pilot symbols spread with a known CDMA spreading 
code. 
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5. The Doppler spread estimator of claim 1 further comprising third 
processing means operatively associated with the receive means for stationarizing the digital 
samples prior to estimating spectral densities of the digital samples. 

6. The Doppler spread estimator of claim 5 wherein said received signal 
5 represents a sequence of transmitted symbols and said third processing means comprises a 

multiplier for mixing the digital samples with a complex conjugate of the transmitted symbols. 

7. The Doppler spread estimator of claim 6 wherein said complex 
conjugate of the transmitted symbols is derived from a symbol estimator. 

8. The Doppler spread estimator of claim 5 wherein said third processing 
10 means further comprises a low pass filter. 

9. The Doppler spread estimator of claim 1 wherein said first processing 
means develops vector data representing spectral components. 

1 0. The Doppler spread estimator of claim 1 wherein said second processing 
means comprises a multichannel correlator. 
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11. The Doppler spread estimator of claim 10 wherein said correlator 
calculates correlations between estimated spectral densities and weighting functions 
corresponding to the hypotheses of the plural Doppler spread values. 

12. The Doppler spread estimator of claim 11 wherein the second 
5 processing means includes a selector for selecting the hypothesis producing a highest 

correlation value. 
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13. A receiving apparatus Doppler spread estimator comprising: 

a receiver receiving signals over a transmission channel that varies over time 
due to fading and converting the signals to digital samples for processing; 

an estimator operati vely associated with the stationarizer for estimating spectral 
5 densities of the stationarized digital samples; and 

a correlator operatively coupled to the estimator for correlating the spectral 
density estimates with functions corresponding to hypotheses of plural Doppler spread values 
to select one of said Doppler spread values as an estimate of actual Doppler spread. 

4 14. The receiving apparatus Doppler spread estimator of claim 13 wherein 
10 said receiver is adapted to receive Code-Division Multiple Access (CDMA) signals. 

15. The Doppler spread estimator of claim 14 wherein channel coefficients 
are directly obtained from the CDMA signals. 

16. The receiving apparatus Doppler spread estimator of claim 13 wherein 
said receiver is adapted to receive signals including pilot symbols spread with a known CDMA 

1 5 spreading code. 

17. The Doppler spread estimator of claim 13 further comprising a 
stationarizer operatively associated with the receiver for stationarizing the digital samples prior 
to estimating spectral densities of the digital samples. 
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1 8. The receiving apparatus Doppler spread estimator of claim 1 7 wherein 
said received signal represents a sequence of transmitted symbols and said stationarizer 
comprises a multiplier for mixing the digital samples with a complex conjugate of the 
transmitted symbols. 

5 19. The receiving apparatus Doppler spread estimator of claim 1 8 wherein 

said complex conjugate of the transmitted symbols is derived from a symbol estimator. 

20. The receiving apparatus Doppler spread estimator of claim 1 8 wherein 
said stationarizer further comprises a low pass filter. 

2 1 The receiving apparatus Doppler spread estimator of claim 1 3 wherein 
10 said estimator develops vector data representing spectral components. 

22. The receiving apparatus Doppler spread estimator of claim 1 3 wherein 
said correlator comprises a multichannel correlator. 

23 The receiving apparatus Doppler spread estimator of claim 22 wherein 
said correlator calculates correlations between estimated spectral densities and weighting 
15 functions corresponding to the hypotheses of the plural Doppler spread values. 
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24. The receiving apparatus Doppler spread estimator of claim 23 wherein 
the correlator includes a selector for selecting the hypothesis producing a highest correlation 
value. 
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25. The method of estimating Doppler spread of a received signal, 

comprising the steps of: 

receiving the signals over a transmission channel that varies over time due to 

fading and converting the signals to digital samples for processing; 
5 estimating spectral densities of the stationarized digital samples; and 

correlating the spectral density estimates with functions corresponding to 
hypotheses of plural Doppler spread values to select one of said Doppler spread values as an 
estimate of actual Doppler spread. 

26. The method of claim 25 wherein said receiving step receives Code 
10 Division Multiple Access signals. 

27. The Doppler spread estimator of claim 26 wherein channel coefficients 
are directly obtained from the CDMA signals. 

28. The method of claim 25 wherein said receiving step receives signals 
including pilot symbols spread with a known CDMA spreading code. 

j 5 29. The Doppler spread estimator of claim 25 further comprising third 

processing means operatively associated with the receive means for stationarizing the digital 
samples prior to estimating spectral densities of the digital samples. 
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30. The method of claim 25 wherein said received signal represents a 
sequence of transmitted symbols and said stationarizing step comprises multiplying the digital 
samples with a complex conjugate of the transmitted symbols. 

31. The method of claim 30 wherein said complex conjugate of the 
5 transmitted symbols is derived from a symbol estimator. 

32. The method of claim 30 wherein said stationarizing step comprises low 
pass filtering the mixed samples. 

33. The method of claim 25 wherein said estimating step develops vector 
data representing spectral components. 

10 34. The method of claim 25 wherein said correlating step comprises using 

a multichannel correlator. 

35.- The method of claim 34 wherein said correlator calculates correlations 
between estimated spectral densities and weighting functions corresponding to the hypotheses 
of the plural Doppler spread values. 

15 36. The method of claim 3 5 wherein the correlating step includes selecting 

the hypothesis producing a highest correlation value. 
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